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We demonstrate a one-dimensional optical lattice clock with a spin-polarized fermionic 
isotope designed to realize a collision-shift-free atomic clock with neutral atom ensembles. 
To reduce systematic uncertainties, we developed both Zeeman shift and vector light-shift 
cancellation techniques. By introducing both an H-maser and a Global Positioning System 
(GPS) carrier phase link, the absolute frequency of the 1 S (F = 9/2) - 3 P (F = 9/2) clock 
transition of the 87 Sr optical lattice clock is determined as 429,228,004,229,875(4) Hz, where 
the uncertainty is mainly limited by that of the frequency link. The result indicates that the 
Sr lattice clock will play an important role in the scope of the redefinition of the "second" 
by optical frequency standards. 
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1. Introduction 

The rapid development of research on optical frequency measurement based on femtosec- 
ond combs 1 -*' 2 ) has stimulated the field of frequency metrology, especially research on high- 
performance optical frequency standards. Optical frequency standards based on single trapped 
ions 3 ) -6 ) and ultracold neutral atoms in free fall 7 )' 8 ) have provided record levels of performance 
that approach those of the best Cs fountain clocks with a fractional frequency uncertainty of 
below 1 x 10~ 15 . 9 ) 

We have proposed a novel approach named an "optical lattice clock", in which atoms 
trapped in an optical lattice potential serve as quantum references. 10 )~ 12 ) The sub- wavelength 
localization of a single atom in each lattice site suppresses the first-order Doppler-shift and 
collisional- frequency-shift while it provides a long interrogation time of over 1 s. The light shift 
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induced by the trapping field can be precisely canceled out by carefully tuning the lattice laser 
wavelength to the "magic wavelength". 10 )' n )> 13 )~ 15 ) 

Recently, a higher order light shift, which is not canceled out at the "magic wavelength" 
and that imposes an uncertainty limit on the lattice clock scheme, was observed in a Sr lattice 
clock and revealed to affect only the accuracy at a level below 10~ 18 . 16 ) Optical lattice clocks 
have already achieved a linewidth that is one order of magnitude narrower 12 )' 17 ) than that 
observed for conventional neutral-atom optical clocks. 7 )' 8 ) The very high potential stability 
as well as accuracy of this scheme would permit the measurement of a fractional frequency 
difference at the 1CP 18 level in a minute, 11 ) which may open up new applications for ultra- 
precise metrology, such as the search for the time variation of fundamental constants, 4 ) the 
real time monitoring of the gravitational frequency shift and the redefinition of the "second" . 

Until now, optical lattice clocks have only been demonstrated with one-dimensional (ID) 
optical lattices employing spin-unpolarized fermions 12 )' 15 )' 18 )' 19 ) or bosons. 17 ) While these 
experiments have clearly demonstrated the advantage of Lamb-Dicke confinement 20 ) provided 
by an optical lattice, collisional-frequency shifts should exist because of the relatively high 
atomic densities of up to 10 11 /cm 3 trapped in a single lattice site. 15 )' 21 ) This collision shift 
would ultimately be a fatal accuracy problem for ID optical lattice clocks as witnessed for 
such neutral atom based clocks as Cs fountain clocks 9 ) or Ca optical clocks, 7 )' 8 ) in which 
collision shifts dominate their uncertainty budgets. Even in the presence of other particles, 
it has been predicted 22 ) and demonstrated 23 )' 24 ) in the RF transition that the collisional 
frequency shifts can be suppressed through the Pauli exclusion principle 25 ) -27 ) by employing 
ultracold spin-polarized fermionic atoms. 

In this paper, we demonstrate, for the first time, a ID optical lattice clock with ultracold 
spin-polarized fermionic atoms, which, in principle, would realize collisional-shift-free atomic 
clocks. In addition, the Zeeman shift and the vector light shift cancellation technique have 
been introduced to further improve the clock accuracy. Furthermore, an improved frequency 
measurement based on an H-maser and a Global Positioning System (GPS) carrier phase link 
is discussed in detail. 

The absolute frequency of the transition for the Sr lattice clock was first determined to 
be 429,228,004,229,952(15) Hz using a Cs clock referenced to the SI second. 12 )' 28 ) Later the 
JILA group measured the frequency and found it to be 429,228,004,229,869(19) Hz. 18 ) The 
measurement results obtained by the two groups were in poor agreement at a level of three 
times the combined uncertainties. 

In order to resolve this inconsistency, we have improved the absolute frequency measure- 
ments based on an H-maser linked to UTC (NMIJ) using GPS carrier phase signals. The UTC 
(NMIJ) is in turn linked to international atomic time (TAI). The Allan standard deviation 
is obtained for the Sr lattice clock and is found to reach 2 x 10~ 15 at an averaging time of 
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Fig. 1. Energy levels for Sr atoms. Spin-polarized ultracold Sr atoms were prepared by using the 
1 S'o — 1 P\ transition at A = 461 nm and the 1 Sq — 3 Pi transition at A = 689 nm. The first-order 
Zeeman shift on the clock transition at A = 698 nm was eliminated by averaging the transition 
frequencies f±g/2, corresponding to the 1 So(F = 9/2, mp = ±9/2) — 3 Po{F = 9/2, mp = ±9/2) 
clock transitions, respectively 

1300 s. The newly obtained absolute frequency in this work is 429,228,004,229,875 Hz, with 
an uncertainty of 4 Hz. This frequency value differs from that of our previous measurement 
by five times the combined uncertainty but falls within the uncertainty of the JILA value. 
We reported the preliminary results of our improved frequency measurement at CLEO/QELS 
20 06. 29 ) Later we learned that the SYRTE group had also reported a measured frequency value 
for the Sr lattice clock of 429,228,004,229,879(5) Hz on the arXiv 19 ) during the CLEO/QELS 
conference. There is good agreement between the measurement results obtained by the three 
groups. 

2. Method 

2.1 Experimental setup 

Figure 1 shows relevant energy levels for 87 Sr atoms. Strontium atomic beams effused from 
an oven heated at 800 K were decelerated and magneto-optically trapped on the 1 Sq — 1 P\ 
transition at A = 461 nm down to a few mK. The atoms were further cooled and trapped using 
the Dynamic Magneto-Optical Trapping (DMOT) scheme 21 ) on the 1 Sq — 3 Pi transition at A = 
689 nm. Ultracold Sr atoms at a few fiK were then loaded into a ID optical lattice consisting 
of a standing wave laser operated at the "magic wavelength" of Al = 813.428 nm. 12 )' 16 ) The 
lattice laser was focused into an e -2 beam radius of 30 /xm. At the anti-node of the standing 
wave a peak power density of /l =10 kW/cm 2 was obtained, which gave axial and radial trap 
frequencies of v x = 40 kHz and v y (K> v z ) = 350 Hz, respectively. 
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Fig. 2. Schematic diagram of the experimental setup. AOM, acousto-optic modulator; PC, personal 
computer; DBM, double balanced mixer; ULE, ultra-low expansion; GPS, Global Positioning 
System; UTC, coordinated universal time; NMIJ, National Metrology Institute of Japan; TAI, 
international atomic time. A beat signal was used to stabilize the n-th comb component to the 
Sr transition frequency. 



Figure 2 shows a schematic of the experiment. A clock laser operating at Ao = 698 nm 
was frequency-stabilized to a high-finesse ULE cavity with its finesse of 430,000 to reduce 
the laser linewidth to about 20 Hz. The clock laser was steered into a Sr optical lattice 15 ^' 30 ) 
through a polarization-maintaining single-mode optical fiber, where a fiber noise canceller 31 ) 
was installed. An acousto-optic modulator (AOM) operating at about 40 MHz was used for 
both frequency control and intensity switching to produce an excitation 7r-pulse. The beam 
waist of the clock laser was 300 /jm to reduce intensity inhomogeneity. 
2.2 Spectroscopy of trapped atoms 
2.2.1 Spin polarization and cold collision suppression 

In ID optical lattice clocks, collisional frequency shifts should occur in each lattice site, 
where typically a few tens of atoms were trapped, corresponding to an atomic density of 
10 n /cm 3 . Its suppression, therefore, is an important matter to be resolved if we are to realize 
accurate atomic clocks. 

In ultracold collisions, in general, the partial waves that contribute to the collisions are 
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Fig. 3. Atoms loaded into a ID optical lattice were spin-polarized by a er^ polarized optical-pumping 
laser operating on the 1 Sq(F = 9/2) — 3 Pi(F = 9/2) transition at 689 nm in the presence of a 
bias magnetic field of B . A clock laser with a wave- vector k was introduced along the ID lattice 
axis and excited the 1 Sq(F — 9/2) — 3 Pq(F — 9/2) clock transition. The electric field vectors for 
both the lattice laser El and the clock laser Eq were parallel to the bias magnetic field Bq. 



greatly limited by the centrifugal barrier for the relative angular momentum I of collision 
pairs. The effective potential is given by, 



C 6 h 2 1(1 + 1) 



(1) 



where r is the inter-atomic distance, Cq the van der Waals constant, h the Planck constant, and 
[i the reduced mass of collision pairs. For example, the p-wave (I = 1) barrier for Sr atoms in the 
1 So ground state is estimated to be about 96 fiK assuming Cq = 3103 a.u. (atomic units), 32 ) 
which is well above the atomic temperature of 3 used in this experiment. In addition, 
even partial waves, such as the s-wave, are not allowed for spin-polarized fermions due to the 
anti-symmetrization of the wavefunction. Therefore we expect the collisional frequency shift 
to be effectively suppressed by spin-polarizing ultracold fermionic 87 Sr atoms. Furthermore 
it has been demonstrated 23 )' 24 ) in an RF spectroscopy of ultracold fermions that the atoms 
remain identical and cannot interact in the s-wave regime in the coherent transfer process. 

Spin polarization was performed as follows: In a presence of a bias magnetic field |Bo| = 
50 mG applied perpendicular to the lattice beam axis as shown in Fig. 3, a circularly 
polarized pumping-laser resonant to the 1 So(F = 9/2) — 3 P±(F = 9/2) transition was applied 
for 50 ms to optically-pump the atomic population to the mp = ±9/2 Zeeman substates 
in the 1 Sq(F = 9/2) ground state, respectively. With this optical pumping scheme, heating 
of spin-polarized atoms can be minimized, as the mp = ±9/2 final states are the "dark 
states" with this pumping laser. The typical atomic temperature was about 3 //K after spin- 
polarization, and more than 95 % of the atoms were transferred to the stretched mp = ±9/2 
state depending on the helicity of the pumping laser. 

Then the bias magnetic field was increased to |Bo| = 1.78 G to completely resolve adjacent 
Zeeman components in the 1 S (F = 9/2, m F = ±9/2) - 3 P (F = 9/2, m F = ±9/2) clock vr- 
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transition (see Fig. 1). This clock laser with its wave-vector k(= (k x ,k y ,k z )) was introduced 
along the lattice beam. The electric field vector of the clock laser Ec and the lattice laser El 
were parallel to the bias magnetic field Bo as shown in Fig. 3. 

2.2.2 Dephasing of Rabi oscillation 

The Rabi frequency fi n , n for a trapped atom in the (n x ,n y , n z )-th vibrational level in the 
lattice potential is described as, 33 ) 

tt n ,n = X| \(nj\e^p(ikjj)\nj}\, (2) 

j=x,y,z 

where is the Rabi frequency of the electronic transition (clock transition), \rij) the wave 
function of the harmonic oscillator state of an atom in the lattice potential along j = x,y,z 
direction, and (x, y, z) the displacement of the atom. Here we assumed an elastic component 
(Arij = 0) in the vibrational transition that is relevant to the clock signal. Using Laguerre 
polynomial L n (x), each matrix element, say j = x component, can be written as, 33 ) 

(n x \ exjp(ik x x)\n x ) = exp(-7^/2)L n ^ {rfc), (3) 

where ij x = k x XQ is the Lamb-Dicke parameter, xq = a/ h / (2mu x ) / (2tt) the spatial extent of 
the ground state wavefunction along x-axis, and m the mass of an trapped atom. The Lamb- 
Dicke parameter for x-direction was r) x = 0.34. Similarly we estimated that for the radial 
direction rj y = k y yo and r\ z = k z zo to be less than 0.05, where we assumed the uncertainty in 
aligning the clock laser with respect to the lattice laser to be less than 12 mrad. 

These finite Lamb-Dicke parameters rjj gave rise to the variation of Rabi frequencies 
0„ in that depended on the vibrational states \rij) of the atoms through Laguerre polynomial 
L nj (rfj). The thermal occupation of the vibrational states was typically n x ~ 1.1 and fi y w 
fi z ~ 180 for the trapped atom temperature. The thermal distribution of the vibrational states 
introduced a vibrational-state-dependent Rabi frequencies for each trapped atom, leading 
to the dephasing of the Rabi oscillation as shown in Fig. 4(a). As a result, the excitation 
probability for the 2-ms-long ir clock pulse was degraded to about £ « 0.8, as deduced by a 
fit in Fig. 4(a). In this measurement, the clock laser was tuned to the atomic resonance with 
an intensity of 25 mW/cm 2 . 

This dephasing of the Rabi oscillations can be suppressed by applying sideband cooling 
to the axial motion and Doppler cooling to the radial motion, 21 ) and by improving the beam 
overlap between the clock and lattice lasers. 

2.2.3 Normalization of the spectrum 

The excitation of the clock transition was observed by measuring the laser induced fluores- 
cence by driving the 1 S$ — 1 P\ cyclic transition as shown in Fig. 1. This fluorescence intensity 
Is was proportional to the atom number Ns remaining unexcited in the 1 Sq ground state. In 
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Fig. 4. (a) Rabi oscillation of atoms in the ID optical lattice. The excited atom fraction n was 
measured as a function of the duration of the clock laser resonant to the 1 So(F = 9/2) — 3 Pq(F = 
9/2) transition, (b) Typical clock spectrum obtained with a 10-ms-long 7r-pulse clock laser. The 
excited atom fraction k was measured as a function of the detuning of the clock laser. 



order to normalize the excited atom fraction, we measured the number of atoms Np in the 
3 Po state. After blowing out the unexcited atoms in the ground state by irradiating the laser 
resonant with the 1 So — 1 P± transition, the atoms in the 3 Po state were deexcited to the 1 So 
ground state by the 0.5-ms-long clock laser with a pulse area of it. Then we measured the laser 
induced fluorescence intensity Ip on the 1 Sq — l P\ transition. This fluorescence intensity was 
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proportional to £Np, where £ ~ 0.8 was the (de)excitation efficiency for the ir clock-pulse 
that was given by the inhomogeneous Rabi frequencies as discussed previously. The excited 
atom fraction k was thus calculated as, 

K= N ? = !z ( 4) 

Ns + Np ZIs + Ip' K> 
Figure 4(b) shows the typical clock transition excited with the 7r-pulse clock laser with an 

intensity of 1 mW/cm 2 . The excited atom fraction n = k(5v) was measured as a function of 

the clock laser detuning 5v from the atomic resonance. The duration of the clock excitation 

pulse was 10 ms, and a nearly Fourier limited linewidth of 80 Hz was observed. 

2.2.4 Frequency stabilization and Zeeman shift cancellation 

Frequency stabilization 34 ) of the clock laser to the spectrum center was realized by feedback 
control of the AOM frequency (see Fig. 2) using the error signal f e rr(t n ) obtained by a digital 
servo loop as, 

jcrr(Wl) = /err(*n) + Sf(t n ). (5) 

Here 5f(t n ) is the correction signal measured in the n-th interrogation period at t = t n as, 

*/ = 7X * (+ ^V 7/2) - W 
where 7 ~ 80 Hz is the full width at half maximum (FWHM) linewidth of the observed 
spectrum, and re(±7/2) is the atom excitation probability near the side slopes of the Rabi 
excitation spectrum with respect to the stabilized line center at t = t n . In this servo loop, 
the drift rate of the clock laser frequency was evaluated in advance and was fed forward 
to minimize servo errors occurring during the frequency stabilization. We estimate that the 
servo-error would be of the order of 0.1 Hz. The application of a second order integrator, in 
addition to the currently used first order integrator, will further reduce servo errors. 

The measured spectrum was Zeeman shifted by a bias magnetic field of Bo- As depicted 
in Fig. 1, the Zeeman splitting is smaller in the s Pq(F = 9/2) state than in the 1 Sq(F = 9/2) 
state due to hyperfine mixing in the 3 Pq(F = 9/2) excited state. This introduces a first order 
Zeeman shift of nip x 106 Hz/G 11 ) for the 7r-transition excited from the mp sublevel in the 
1 So(F = 9/2) ground state. This linear Zeeman shift can be eliminated by averaging the 
two transition frequencies f±g/2, corresponding to the 1 Sq(F = 9/2, mp = ±9/2) — 3 Pq(F = 
9/2, mp = ±9/2) transition frequencies, respectively. The transition frequency free from the 
first order Zeeman shift is thus given by, 

, Z+9/2 + 7-9/2 . . 

Jo = 2 • 

Since a single measurement took 1 s or less, a cycle time t c = t n — t n -\ of nearly 2 s 
was required to determine one of the Zeeman components f±g/2- The cycle time required for 
cooling, trapping, and interrogating atoms in the lattice, was not optimized in this experiment. 
A further reduction of the cycle time to less than 1 s will be feasible in a future experiment. 
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2.3 Frequency measurement 

To measure the absolute frequency of the Sr lattice clock, a frequency comb system was 
brought from NMIJ, AIST in Tsukuba to the University of Tokyo, where the optical lattice 
clock was operated. Instead of the commercial Cs clock that was used as a local time base in 
our previous frequency measurement, 12 **' 28 -* we moved an H-maser (Kvarz, Model CH1-75A) 
from Tsukuba to Tokyo to reduce the measurement uncertainty. 

In the frequency comb, the frequency of the n-th comb component is expressed as 

fn = nX / rC p + /cEO, (8) 

where / rep is the repetition rate of the laser pulse and /ceo is the carrier-envelope offset 
frequency. 2 ^ When / rep and /ceo are precisely controlled, the comb works as a "frequency 
linker" that connects optical and radio frequencies. The control ports of / rep and /ceo i n the 
Ti:sapphire laser are the lengths of the laser cavity and the pump laser power, respectively. 
Our comb system is described in detail elsewhere. 28 )> 35 ) 

In this experimental configuration, since the laser light after the AOM was a pulsed-light 
used for the Sr spectroscopy, its frequency could not be directly measured with the frequency 
comb. The frequency relations in the frequency measurement are shown in Fig. 2. We first 
measured the beat frequency /b = \f c — f n \ between the clock laser f c and the n-th tooth 
of the comb /„ with a photo-diode. We then electronically mixed the beat note /b with 
/aom = /err + 40 MHz with a double balanced mixer (DBM) and extracted the frequency 
component, 

fb = |/c + /AOM - fn\, (9) 

which corresponded to the + 1 order light diffracted by the AOM. This frequency is equal to a 
beat frequency between the Sr-transition frequency /s r = / c + /aom and the n-th tooth of the 
comb f n . In our measurement scheme, / b was used to phase- lock the n-th comb component 
to the Sr clock transition by feedback controlling the cavity length of the mode-locked laser. 

In this way, the whole comb was locked to the Sr clock transition, which means that the 
stability of each comb component and / rep follows that of the Sr clock transition. / rep was 
measured against the H-maser as follows. / rep was observed at about 793 MHz and down 
converted to a frequency of / r ' cp = 800 MHz — / rep ~7.3 MHz by using a DBM and a low-pass 
filter. The frequency of 800 MHz was generated by using a synthesizer with a fixed frequency 
of exact 800 MHz. /' was measured and recorded with a universal counter (Agilent, Model 
53132A). All the synthesizers and counters used in this experiment were phase locked to the 
H-maser through a distribution amplifier. Finally, /g r was calculated by using the equation: 

fsr = n x (800MHz - / r ' ep ) + / CEO ± K, (10) 

where / r ' ep was measured by the counter, while /ceo and / b were set at fixed frequencies using 
digital phase-lock loops. The integer n was simply determined by solving eq. (8) for n and 
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requiring n to be an integer, since we already know /g r with an uncertainty much smaller 
than / rep . The sign of /■£ in eq. (10) was determined by changing / rep slightly and observing 
the variation of the // in the experiment. 

To calibrate the frequency of the H-maser, GPS carrier phase receivers (Javad, Model 
$Texon-GGD) were employed at both sites, the University of Tokyo and NMIJ. In our previ- 
ous measurement, 12 )> 28 ) a GPS disciplined oscillator was introduced to link the local Cs clock 
to the GPS time. 36 ) To further improve the link precision, in the present experiment, the 
H-maser was calibrated based on UTC (NMIJ) using the GPS carrier-phase technique with 
the analysis software "GIPSY". 37 ) The distance between Tokyo and Tsukuba is about 51 km. 
The relationship between the UTC (NMIJ) and TAI can be found in the monthly reports of 
the Circular-T of the Bureau International des Poids et Mesures (BIPM). 38 ) 

3. Experimental results 

3. 1 Stability evaluation of the lattice clock locked to the two Zeeman components 

In order to evaluate the clock stability at optical frequencies, the optical lattice clock 
was alternately stabilized to two Zeeman components, i.e., the 1 Sq(F = 9/2, m F = ±9/2) — 
3 Pq(F = 9/2, vriF = ±9/2) clock transition (see Fig. 1), corresponding to the transition 
frequencies f±g/2, respectively. The duration of the clock excitation pulse was 10 ms, and a 
nearly Fourier limited linewidth of 80 Hz was observed. Four successive measurements were 
used to lock the clock laser frequency to the /+g/2 and /_g/2 transition frequencies. 

Error signals / crr that were fed back to the AOM (see Fig. 2) are shown in Fig. 5(a) as a 
function of elapsed time. The top and bottom curves correspond to the case when the clock 
laser frequency is locked to the j+9/2 (top) and f-g/2 (bottom), respectively. The curve in the 
middle shows the average of these two Zeeman components fo (see eq. (7)), which provided 
the first-order Zeeman-shift-free transition frequency as discussed previously. The slope of 
these curves compensated the ULE cavity drift rate, which was found to be —0.16 Hz/s. 

The two clock frequencies f+9/2 and f-g/2 can De regarded as output signals generated 
by two independent optical clocks locked to different Zeeman components. The optical beat 
note of these two clocks A/ = j+9/2 — f-9/2 can be- evaluated by the offset of these two error 
signals, namely the top and bottom curves in Fig. 5(a), and was used to evaluate the stability 
of lattice clocks. Figure 5(b) shows the Allan standard deviation, which was measured for 
up to 7600 s. Until 10 s, when the servo loop started to work, the Allan standard deviation 
increased. After that the deviation started to decrease with cr(r) = 8 x lO -14 /^. The floor 
of the Allan deviation was not observed during this measurement for a time of over 2 h. 

3.2 Frequency measurement and stability evaluation 

Figure 6 shows the measured frequencies of the Sr clock transition against the H-maser. In 
Fig. 6(a), the clock laser was stabilized to the ^oCF = 9/2, m F = -9/2)- 3 P (^ = 9/2, m F = 
—9/2) Zeeman component with f-g/2, while in Fig. 6(b) the clock laser was stabilized to the 
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Fig. 5. (a) Error signals / C1T fed back to the clock laser stabilized to the ULE cavity as a function 
of time. The three curves show the error signals corresponding to f+g/2 (top), fo (middle), and 
f-9/2 (bottom), respectively, fo provided the first-order Zeeman-shift-free clock transition fre- 
quency The slope of these curves indicated the ULE cavity drift rate of —0.16 Hz/s. (b) Allan 
standard deviation evaluated by the optical lattice clock independently locked to the two Zeeman 
components as shown in (a). 



center of the clock transition /o by using the spectra of both Zeeman components, i.e., f+9/2 
and /-_9/2- The gate time of the universal counter was set at 10 s. The scatter of the measured 
frequencies in both figures mainly results from the instability of the H-maser. 

The measurement time for Figs. 6(a) and (b) were 9140 s and 7610 s with average frequen- 
cies of 429,228,004,229,019.5 and 429,228,004,229,867.0 Hz, respectively. Note that systematic 
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Fig. 6. Variation in the Sr-stabilized laser frequency measured against the H-maser as a function 
of time. The averaging time was 10 s. Systematic corrections are not applied, (a) Stabilized to 
the 1 S {F = 9/2, m F = -9/2) - 3 P {F = 9/2, m F = -9/2) Zccman components of /_ 9/2 ; (b) 
stabilized to the line center /q. 



corrections are not applied in these figures. A total of ten measurements were performed on 
different days over a period of more than one month. The measurements in Figs. 6(a) and (b) 
are denoted as measurement #2 and #10, respectively. 

The Allan standard deviations were calculated and indicated in Fig. 7: A solid curve with 
filled triangles and a dashed curve with filled circles correspond to the measured frequencies 
indicated in Figs. 6(a) and (b), respectively. Both curves show similar trend for the Allan 
deviation of the Sr-stabilized laser as shown in Fig. 5(b). For comparison, the Allan standard 
deviation of the H-maser is also shown in Fig. 7 as a solid curve. We notice that, for the short 
term (averaging time between 1 and 20 s), the measured stability was basically limited by the 
H-maser. For the long term (averaging time longer than 20 s), the measured Allan deviation 
shows the fractional frequency instability of the Sr lattice clock, which is consistent with the 
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Fig. 7. Allan standard deviation of the Sr lattice clock measured based on an H-maser. A solid curve 
with filled triangles and a dashed curve with filled circles correspond to the measured frequencies 
indicated in Figs. 6(a) and (b), respectively. A solid curve below indicates the Allan standard 
deviation of the H-maser. 

Allan standard deviation derived from the optical beat note shown in Fig. 5(b). The Allan 
standard deviation of the Sr lattice clock has reached a level of 2 x 10 -15 at an averaging time 
of 1300 s. This means that the statistical error of each measurement is reduced to about 0.9 
Hz after a 1300 s averaging time. 

4. Frequency corrections and uncertainties 

4-1 Frequency link 

To reach the final measurement result, systematic corrections were applied to the obtained 
averaged frequencies. In the present measurement, the main correction is from the H-maser, 
which was calibrated based on the UTC (NMIJ) using the GPS carrier-phase technique. 

Figure 8(a) shows a typical recording of the time difference between the UTC (NMIJ) re- 
alized in Tsukuba and the H-maser located at the University of Tokyo. The data was recorded 
every 300 s over 2 days. The slope of the curve indicates the fractional frequency offset between 
the UTC (NMIJ) and the H-maser, which was calculated to be —1.33 x 10~ 14 using a linear 
fitting. The minute variations in the measured data around the fitted solid line are considered 
to be mainly caused by the atmospheric fluctuation along the link path. Figure 8(b) shows 
the Allan standard deviation of the GPS link calculated from the measured data shown in 
Fig. 8(a). The Allan standard deviation was 6.6 x 10 -14 for a 300 s averaging time, which 
improved to 3 x 10~ 15 after a 30000 s averaging time. 
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10 3 10 4 

(b) Averaging Time (s) 

Fig. 8. (a) Time difference between the UTC (NMIJ) and the H-maser measured by the GPS link 
between Tokyo and Tsukuba using carrier-phase signals. The slope indicates the fractional fre- 
quency difference between the UTC (NMIJ) and the H-maser. (b) The frequency stability of the 
GPS link. 

During these two days, we performed measurements #7 and #8. Based on the results 
of the GPS link, the correction from the H-maser was calculated to be +5.7 Hz as listed in 
Table I. The uncertainty of the correction is basically limited by the stability of the GPS link 
obtained within the measurement time (also listed in Table I). 

The uncertainties of the H-maser correction were calculated to be 3.9 and 3.4 Hz for mea- 
surements #7 and #8, respectively. The corrections and uncertainties of all ten measurements 
contributed from the calibration of the H-maser are listed in Table I. 

Note that the correction from the H-maser gradually varied from —10.8 to + 10.7 Hz during 
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Tabic I. Corrections and uncertainties of the H-maser calculated based on the GPS link. 



Meas. # 


Measurement 


Correction 


Uncertainty 




time (s) 


(Hz) 


(Hz) 


1 


910 


- 10.8 


10.7 


2 


9140 


- 10.8 


4.7 


3 


4020 


- 10.8 


6.4 


4 


3460 


+ 2.7 


5.1 


5 


4530 


+ 3.8 


5.1 


6 


6580 


+ 3.8 


4.7 


7 


5480 


+ 5.7 


3.9 


8 


6470 


+ 5.7 


3.4 


9 


9320 


+ 10.7 


3.0 


10 


7610 


+ 10.7 


3.4 



the 38-day measurement period, which was due to a fractional frequency drift of 6.6 x 10 for 
the H-maser. The fractional frequency drift of the H-maser at the NMIJ was about 1 x 10~ 14 
over one month. The increased frequency drift at the University of Tokyo is considered to be 
due to the short setting time and the difference in the laboratory environment after moving 
the H-maser. 

The results of the frequency link between the UTC (NMIJ) and the TAI can be found in 
Circular T, which is published on the BIPM web page. 38 ) From Circular T 221 and 222, the 
correction and the uncertainty were calculated to be and 1.4 Hz, respectively (see Table II). 
UTC (NMIJ) has been very stable during the past several months. 

A frequency counter usually counts a frequency with systematic and statistical errors. The 
universal counter (Agilent 53132A) used in the present experiment was measured and found 
to have a fractional systematic error (correction) of +7.4 x 10~ 13 and a fractional statistical 
error (uncertainty) of 2.0 x 10~ 13 . Since the counter is used to read the / r ' ep at 7.3 MHz instead 
of the / rep at 793 MHz, the correction and the uncertainty applied to the / rcp is —6.8 x 10~ 15 
and 1.8 x 10~ 15 , which corresponds to a correction of —2.9 Hz and an uncertainty of 0.8 Hz in 
the Sr optical frequency (see Table II). The systematic correction and statistical uncertainty 
of the counter can be further reduced by increasing the gate time in the measurement. By 
observing the higher order of / rep with a fast detector and increasing the frequency of the 
synthesizer, the effect of the counter errors in the frequency measurement can also be further 
reduced. 
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Table II. Systematic corrections and uncertainties for the Sr lattice clock. The other corrections or 
uncertainties that are less than 0.1 Hz are omitted from the list. Corrections for the GPS link 
were applied to each data set as summarized in Table I. 



Effect 


Correction 


Uncertainty 




(Hz) 


(Hz) 


[Frequency Measurement] 






r~\ OC V 1 TT^ 1 TTTi/~l/"lvnv /TT TM 

GFb ImJqlokyo - U 1L(JNM1J)J 




3.8 


T~\ ... l- l Tt Tmn /ATA /TT T \ rp A t"| 

hreq. lmk[U IC+lNMUj-lAl] 





1.4 


Counter reading 


— z.y 


0.8 


[Sr Optical Lattice Clock] 






Blackbody radiation shift 


+2.4 


0.2 


Gravitational shift 


-0.5 


0.1 


2nd order Zeeman shift 


+0.736 


< io- 3 


Scalar light shift 





0.5 


Probe laser light shift 


+0.03 


< 10- 3 


Total 


-0.2 


4.2 



4-2 A ID optical lattice clock 

The corrections and uncertainties for the Sr lattice clock are listed in Table II. Relatively 
large corrections come from the blackbody radiation and the gravitational shifts. The black- 
body radiation shift 11 )' 39 ) is calculated to be -2.4 Hz 39 ) for the temperature (T = 301 K) of 
the vacuum chamber in which the optical lattice clock was operated. This temperature was 
slightly higher than room temperature (T = 293 K) because of the heat generated by the 
trapping coils. The uncertainty is estimated assuming the worst case of the chamber temper- 
ature inhomogeneity of AT = ±5 K. For the gravitational shift, the correction is calculated 
using the height of the setup from the geoid surface (11m) while the uncertainty is estimated 
from the uncertainty of the height measurement (2 m). 

As discussed previously, thanks to the application of spin-polarized ultracold fermions, 
the collisional frequency shift can, in principle, be eliminated. However, imperfect spin- 
polarization and inhomogeneous Rabi frequencies £l n ,n, which caused dephasing of Rabi os- 
cillations among trapped atoms, might degrade the Fermi suppression, and introduce a colli- 
sional frequency shift. In spite of this fact, no collisional frequency shift was detected within 
the uncertainty of the present measurements. 

The bias magnetic field |Bo| = 1.78 G applied during the measurements caused a first- 
order Zeeman shift of ±848 Hz for f±g/2, which was eliminated by alternately polarizing atoms 
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into the 1 Sq(F = 9/2, itif = ±9/2) stretched states and averaging two Zeeman components 
/±9/2> as discussed previously. The second order Zeeman shift for this bias magnetic field of 
| Bo | is estimated to be —0.736 Hz. 

It is important to mention that this Zeeman shift cancellation scheme allows us to elim- 
inate the vector light shift that occurs for an elliptically polarized lattice laser, which has 
been regarded as one of the largest uncertainty sources in the error budgets in an optical 
lattice clock. 11 ) In the limit of the negligible mixing of Zeeman sublevels due to the light shift 
perturbation, the vector light shift 40 ) e v is described as the interaction between the effective 
magnetic field <5B, which is proportional to the degree of circular polarization of the lattice 
laser and is parallel to the wave vector of the lattice laser kL , and the magnetic dipole operator 

e v = -/z<5B. (11) 

Therefore, the total magnetic field can be written as B = Bo + 5B. We note that in our present 
configuration <5B is more than thousands times smaller than Bo and is nearly perpendicular 
to Bo- Since the total shift is proportional to rnp|B| and B is assumed to be constant during 
the experiment, the vector light shift differences 5e v for the ±mp Zeeman sublevels can be 
cancelled out when evaluating eq. (7) as in the case of Zeeman shifts. This cancellation allows 
us using the tuf = ±9/2 stretched states in the 1 Sq(F = 9/2) ground state as the clock 
transition as demonstrated in this work, although they have the highest sensitivity for the 
vector light shift. 11 ) 

As regards the tensor light shift, since we interrogated the 1 Sq(F = 9/2, \mp\ = 9/2) — 
3 Pq(F = 9/2, \mp\ = 9/2) transition with a 7r-polarized clock laser, there was no variation in 
the tensor light shift. Moreover, if we define this transition as the "clock transition", there 
should be no correction for the tensor light shift. 

In the present experiment, the scalar light shift was the largest source of uncertainties, and 
it was estimated to be « 0.5 Hz, assuming the magic wavelength of Al = 813.428 nm 16 ) and 
a lattice laser peak intensity of /l = 10 kW/cm 2 . This uncertainty can be further reduced 
to the mHz level 11 -* by observing a narrower clock transition that improves the precision 
with which the "magic wavelength" is determined. With this lattice laser intensity /l, the 
hyperpolarizability effect, which is proportional to I 2 , was at the mHz level. 11 )' 16 ) 
4-3 Absolute frequency of Sr lattice clock 

Figure 9 shows the absolute frequency of the ten measurements after applying the cor- 
rections in Tables I and II. The filled triangle (measurement #1) represents a measurement 
result obtained by measuring frequencies of f+g/2 and f-g/2 separately and by calculating 
an average /o- On the other hand, the filled square (measurement #2) represents a result 
obtained by combining the measured frequency f-g/2 (shown in Fig. 6(a)) and the frequency 
difference Af = f + g/2 — f-g/2 obtained as described previously, i.e., fo = f-g/2 + A//2. The 
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Fig. 9. Absolute frequency measurement of the 1 S (F = 9/2) — 3 Pq(F = 9/2) transition of 87 Sr 
atoms in an optical lattice. We performed 10 measurements over one month. 



filled circles (measurements #3 — #10) represent the results obtained by stabilizing the clock 
laser to the line center of the two Zeeman components /o- 

The error bars in Fig. 9 only include the GPS link uncertainties (listed in Table I) and 
are given as one standard deviation. The weighted average of the ten measurements gives an 
average frequency of 429,228,004,229,875.3 Hz. The standard deviation of the mean is 1.3 Hz. 
However, we would like to point out that we have sometimes observed nonlinear variations in 
the GPS link data, which may cause significant errors in the determination of the correction 
for the H-maser. These errors are included in the corrections for the H-maser but are not easily 
averaged out. Therefore, we took the scatter of the data (the standard deviation) indicated 
in Fig. 9, which is 3.8 Hz, as the uncertainty of the average frequency. 

This uncertainty is then combined with all the other uncertainties listed in Table II and 
the result gives the final combined uncertainty as 4.2 Hz. To the best of our knowledge, this 
is the lowest uncertainty reported for the determination of the absolute frequency of optical 
lattice clocks. 

In Fig. 10, in addition to the present measurement result, we also show the results from 
the JILA 18 ) and the SYRTE 19 ) groups. The weighted average of the frequencies measured by 
the three groups gives an average frequency of 429,228,004,229,876.6 Hz, which is within the 
error bars of the three groups. The standard deviation of the mean is 3.2 Hz. There is good 
agreement between the measurement results of the three groups. This is an important step 
in research on an "optical lattice clock" , and promises to be a significant contribution to the 
field of frequency metrology. 



18/23 



429,228,004. 

^ 429,228,004. 
N 
I 

- 429,228,004. 

o 

g 429,228,004. 

g" 429,228,004. 

^ 429,228,004. 
429,228,004. 



229,900 
229,890 
229,880 
229,870 
229,860 
229,850 
229,840 



J ': A Tokyo-NMIJ SYRTE 
(This work) ^ 



f 



Weighted average = 
429,228,004,229,877 Hz 



Measurement 



Fig. 10. Comparison of the absolute frequencies of the Sr lattice clock measured by three groups. 



5. Discussion 

The absolute frequency of the Sr lattice clock obtained with our improved frequency 
measurement differs from that of our previous measurement by five times the combined un- 
certainty. It is difficult to identify the origin of this difference since the measurement scheme 
has been completely revised. The improved measurement approach has eliminated several 
pitfalls as listed below. 

(1) An H-maser was transported to the University of Tokyo, where the lattice clock was 
operated, for use as a local time base in place of a commercial Cs atomic clock. The 
short-term stability of the frequency measurement was improved by a factor of 50, which 
is extremely useful when evaluating the measurement system. In some cases, the frequency 
instability of the measurement system is related to the systematic errors. Unrevealed 
systematic errors in the previous measurement system under the noise floor caused by the 
instability of the commercial Cs atomic clock may be the main reason for the disagreement. 

(2) A GPS carrier phase link is introduced to calibrate the H-maser through UTC (NMD). 
This link has reduced the uncertainty of the local clock calibration significantly compared 
with the GPS disciplined oscillator used in our previous measurement. 28 ^ 

(3) The clock laser is stabilized to the lattice clock transition in the present experiment. This 
has eliminated the detrimental effects that arise from the timing synchronization and the 
interpolation in the data treatment for spectroscopy and frequency measurement. 

(4) Some elements in the frequency measurement system, including the distribution amplifier, 
have been replaced to reduce measurement instability. 

The whole experiment involves complicated equipment that needs to be transported from 
one place to another. This has prevented the continuous experiments that are needed to trace 
the main improvement of the system step by step. This situation will be improved in future 



19/23 



J. Phys. Soc. Jpn. Full Paper 

by comparing multiple optical lattice clocks at the University of Tokyo or by measuring an 
optical lattice clock at the NMIJ. 

To further reduce the measurement uncertainty, we must reinforce the frequency link be- 
tween the University of Tokyo and the NMIJ in Tsukuba. An optical fiber network that delivers 
optical and microwave frequencies is an effective method for high-precision time transfer. 41 ) 
In terms of a direct link to TAI, a two-way satellite comparison between Japan and Europe 
using the carrier-phase technique is valuable in reducing the uncertainty of the time transfer. 
Of course, a transportable Cs fountain clock should allow us direct access to a highly accurate 
realization of the definition of the second at the local site. 

As for the optical lattice clock, the averaging of the two Zeeman components has signifi- 
cantly reduced systematic uncertainties originating from the Zeeman shift and the vector light 
shift. Although we have employed a spin polarized ultracold fermionic isotope for the clock 
experiment, we found that the inhomogeneous Rabi frequencies among trapped atoms, which 
caused the dephasing of Rabi oscillations, degraded the clock signal intensity and presumably 
the collision suppression due to the Pauli exclusion principle. This problem can be solved by 
reducing the Lamb-Dicke parameters, which can be accomplished by further cooling atoms to 
recoil temperatures or the vibrational ground state, as we have demonstrated. 21 ) Of course, the 
application of 3D optical lattice clocks as in the original proposal 10 - 1 would completely resolve 
these issues related to inhomogeneous Rabi frequencies and collisional frequency shifts. 

6. Conclusion 

We have demonstrated a one-dimensional optical lattice clock with a spin-polarized 
fermionic isotope, which will, in principle, realize a collision-shift-free atomic clock with neu- 
tral atom ensembles. To reduce systematic uncertainties, we have developed Zeeman shift and 
vector light-shift cancellation techniques. By introducing both an H-maser and a Global Posi- 
tioning System (GPS) carrier phase link, the absolute frequency of the 1 Sq(F = 9/2)— 3 Pq(F = 
9/2) clock transition of the 87 Sr lattice clock is determined as 429,228,004,229,875(4) Hz. To 
the best of our knowledge, this provides the most accurate value for the lattice clock. 

In the uncertainty budget of the Sr lattice clock, the scalar light shift is the only term that 
has an uncertainty (fractional 1.2 x 10~ 15 ) larger than that of the Cs fountain clock (fractional 
1 x 10 -15 ). This uncertainty can definitely be reduced by further improving the measurement 
precision of the clock transition. The determination of the absolute frequency of the Sr lattice 
clock at the same level as the state-of-the-art Cs atomic clock should have a great impact on 
discussions related to the redefinition of the "second" . 
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